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We propose a non-perturbative and gauge invariant derivation of the static potential between a 
heavy-quark (Q) and an anti-quark (Q) at finite temperature. This proper potential is defined 
through the spectral function (SPF) of the thermal Wilson loop and can be shown to satisfy the 
Schrodinger equation for the heavy QQ pair in the thermal medium. In general, the proper poten- 
tial has a real and an imaginary part, corresponding to the peak position and width of the SPF. The 
validity of using a Schrodinger equation for heavy QQ can also be checked from the structure of 
the SPF. To test this idea, quenched QCD simulations on anisotropic lattices (a c = 4a T = 0.039fm, 
Na x N z = 20 2 x (96 — 32)) are performed. The real part of the proper potential below the decon- 
finement temperature (7" = 0.78r c ) exhibits the well known Coulombic and confining behavior. At 
(T = 2.33T C ) we find that it coincides with the Debye screened potential obtained from Polyakov- 
line correlations in the color-singlet channel under Coulomb gauge fixing. The physical meaning 
of the spectral structure of the thermal Wilson loop and the use of the maximum entropy method 
(MEM) to extract the real and imaginary part of the proper potential are also discussed. 
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1. Introduction 

Heavy quarkonium at finite temperature is a both intriguing and challenging subject. In par- 
ticular, its in-medium behavior has been proposed to be a prime signal for the creation of the 
quark-gluon plasma (QGP), expected above the critical temperature T c of the deconfinement tran- 
sition [jl|]. Moreover, the advent of RHIC has made it possible to access the QGP in the laboratory 
[Q] and a partial suppression of charmonium has been experimentally confirmed 

In the meantime, theoretical understanding has also progressed significantly. Originally, a non- 
relativistic Schrodinger equation with screened Coulomb potential was considered for the heavy 
QQ pair above T c [jl|]. Later on, more sophisticated potentials were adopted, such as the color- 
singlet free-energy of a static QQ pair and its variants [|J]. A main drawback of these approaches is 
that there is no firm theoretical foundation to use such potentials within the Schrodinger equation 
of a heavy QQ pair in a hot environment. 1 A more direct approach to heavy quarkonium at finite 
temperature T is to extract the spectral function (SPF) from lattice QCD simulations with the help 
of the maximum entropy method (MEM) An unexpected feature found subsequently was that 
charmonium may survive even above the deconfinement transition up to about 2T C [f7], |||. However, 
a transparent understanding of this result has not been obtained so far. 

In light of these circumstances it is imperative to establish a solid connection between the 
Schrodinger approach and the spectral function approach, based on a proper definition of the in- 
medium potential. A possible way to reach this goal was initially proposed in [g]; the proper 
potential to be used in the in-medium Schrodinger equation was defined from the late-time (t) 
behavior of the forward correlator D > (t,R) of a heavy QQ pair separated by distance (/?). It was 
shown by using hard thermal resummation (HTL) techniques at high T that the potential has both 
a real and an imaginary part. The purpose of the present article is to develop this idea further by 
introducing a spectral decomposition of D > (t,R) (or equivalently the thermal Wilson loop) and to 
explore the non-perturbative derivation of the real and imaginary part of the proper QQ potential 
on the basis of lattice QCD simulations. If this program turns out to be successful, a true physical 
understanding of the heavy-quark bound state in the QGP can possibly be obtained. 

2. Formulation 

2.1 Spectral function for a static QQ pair 

We start with the QQ operator M R (t) defined as M R (t) = y{x)TUp{x,y) y(y), where R = |x— yj 
and t = x° = y . The Wilson-line operator Up(x,y) is chosen to connect the points x = (t,x) and 
y = (?,y) by a straight space-like path P. An arbitrary Dirac matrix is denoted by T, e.g. T = 7^ in 
the vector channel and T = iy 5 in the pseudoscalar channel. Let us introduce the forward correlation 
function of the QQ pair, 

D > (t,R) = (M R (t)Ml(0)) (*>0). (2.1) 

In a fully dynamic setting, heavy quarks can propagate in time, hence the relative distance between 
Q and Q will change accordingly. On the other hand, if we consider the infinite mass case (mg — > 

1 This is in contrast to the situation at zero temperature in which a systematic way to derive the potential to be used 
in the Schrodinger equation has been formulated in the framework of non-relativistic QCD [Hi. 
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oo), the spatial separation R reduces to an external parameter. The spectral function associated with 
Eq.([2~ll) reads 

p((0,R) = (1 -e-^)D>((»,R) = l£ \(n\M R (0)\n')\ 2 (e- pE " - e^ E A 8((0 - (E n , - E n )), (2.2) 



where j3 = 1 /T, D > (ft),/?) is the Fourier transform of D > (t,R), and Z is the full partition function 



of the system with gluons, light quarks and heavy quarks. Eq.([2.2|) is antisymmetric in ft) due to 
the bosonic character of the operator Mr. 

We now consider heavy quarks with nearly infinite mass iuq — ► oo, so that they cannot move 
spatially and their relative distance 7? stays fixed. The spectral function can thus be decomposed 
into the sum of the three contributions, p(co,R) = pi(ft),/?) + pn(ft),/?) + pm(co,R) depending on 
the intermediate states given below ("light" implies light quarks and gluons): 





I 


II 


III 


\n>) 


QQ + light 


light 


Q or Q + light 


\n) 


light 


QQ + light 


Q or Q + light 



Since we are interested in the interaction between Q and Q with a spatial separation 7? at the 
same point in time, only the contributions I and II are relevant for our purpose. Moreover, I and II 
are simply related with each other by pu(co,R) = — pi(— ft),/?), due to the anti-symmetric nature 
of the SPF. Without loss of generality, it is therefore sufficient to focus on the contribution I. By 
taking into account the fact that niQ 3> T, we obtain 



p l (a,R) = ^MMRmn')\ 2 8(a-(eAR)-Sn))e- p£ ", 



(2.3) 



where Zo is the partition function without the heavy quarks, \n) and \n') belong to the case I in 
the above Table. Also, we have defined ft) = ft) — 2niQ, £ n '(R) = E n > — 2iuq and e n = E n . Note 
that e n i(R) is the ^-dependent energy of a QQ pair measured from the total rest mass 2mg, while 
e n is /^-independent since no heavy quarks are present in the state \n). Note that they are both 
T-independent by definition. 

The spectral structure of pi(cb,R) contains all the information of the interaction between the 
heavy quarks in the hot medium. To see its connection to the Schrodinger equation for the heavy 
quark system, let us start with the following relation obtained from the definitions of Df" and p\. 



id t Df(t,R) = 2m Q D^(t,R)+e 

If pi (ft), R) has a distinct Breit-Wigner peak at ft) 
reduces to the Schrodinger equation, 



—iZmnt 



e i5 "cdpi((0,R)dC0. 



(2.4) 



(0(R,T) with a half-width £(R,T), Eq.(0) 



id t D>(t,R) 



2m Q + C0(R, T) - iE, (R, r)J Df (t,R) 

2m Q +ReV(R,T)-ilmV{R,T)\Df(t, R) 



(2.5) 



with V(R,T) being the proper heavy-quark potential at finite temperature. On the other hand, if 
there are no well-defined peaks in pi(cb,R), description of the QQ system in terms of the potential 
is not justified. Note that the existence of Im V(R, T) was first pointed out in [Q] in which the late 
time behavior of D > (t,R) was calculated in hard thermal loop resummation at high T. 
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Figure 1: Wilson lines present in the determination of the imaginary-time QQ correlator D(t,R). The loop, 
staple and handles correspond to the case I, II and III defined in Sec. 2. 2, respectively. 



2.2 Thermal Wilson loop and its spectral decomposition 

To evaluate the spectral function defined in Eq.(^3J) in non-perturbative lattice QCD simula- 
tions, we consider the imaginary-time correlator D(t,R) = (T t M#(t)M^(0)) = (M r (t)mI(0)). 

o<T<i/r 

Since the heavy quarks are assumed not to propagate in the spatial directions, their imaginary time 
propagator has the form S#(x — x', t — t') = Se(t — t')<5(x — x') with Se(t — t') being written in 
terms of temporal Wilson-lines. For the full correlator we obtain 



D(t,R) 



-2mQt 



+A_ 



-2m Q (P- 



Tr[*?M)] 
T ' Tr[Wn(z,R)] 



-A ++ e- m ^ Tr[W ma (r,R) + W mb (z,R)} 



(loop) 
(staple) 
(handles) , 



(2.6) 
(2.7) 
(2.8) 



where the numerical coefficients are defined as e.g. A_| = TrpvWrAH] with being the 

projection operator onto the upper and lower components of the Dirac spinor. #[,11,111(1,/?) are the 
Wilson-loops at finite temperature with different topological structure (loop, staple and handle) as 
shown in Fig. 1. 

From their T-dependence and particle content in the intermediate state, it is easy to see that 
loop, staple and handle correspond to the case I, case II and case III respectively. Namely the loop 
contribution Di(z,R) given in Eq.(|2.6|) has the spectral decomposition, 



Di(t,R) = e~ lmQ% \ e- mT p l ((b,R)da). 



(2.9) 

Thus we arrive at the following formula relating the thermal Wilson loop with the spectral function, 

1 



Tr[n(z,R) 



e- aT pi((o,R)d(o. 



(2.10) 



Note that this relation is well-defined in the limit niQ — > 00. 

Lattice QCD simulations of the left hand side of Eq.( 2.10 ) for different values of T and/? allow 
a non-perturbative determination of the spectral function through the inverse Laplace transform, 
e.g. by using the MEM. The peak position and its width can consequently be translated into the 
real and imaginary part of the proper potential. 



3. Quenched QCD results at low and high T 

To study the feasibility of the method proposed above, we perform quenched lattice QCD sim- 
ulations using the plaquette gauge action with j3i at = 7.0 on an anisotropic N% xiV T = 20 3 x (96 — 
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x [fm] R[fm] 

Figure 2: (Left) The thermal Wilson loop W{(z,R) as a function of T at T = 0.78T C . The solid lines are the 
results of the single exponential fit in the intermediate region of T. (Right) Proper potential reconstructed 
from the single exponential fit of the thermal Wilson loop. The gray line here corresponds to a fit based on 
the Cornell-type potential V (R) = c - a/R + OR with a string tension of y/a = 367(5)MeV 



32) lattice. The physical lattice spacing and anisotropy are the same with the first reference in 
[j7|], a a = 4a T = 0.039fm. We adopt the fixed scale method [11] where one varies N z to change the 
temperature. An advantage of this method is that the lattice spacing is the same for all temperatures 
so that the QQ potentials for different temperatures can be directly compared without any adjust- 
ment. The thermal Wilson loop W\(z,R) is calculated as a function of z and 7?. We report here only 
the results at the lowest temperature (T = 0.1ST C ) and the highest temperature (T = 2.33T C ) using 
(125,980) gauge configurations respectively. 

Shown in the left panel of Fig.2 is the on-axis thermal Wilson loop W\ as a function of z for 
different values of R in the low-temperature confinement phase (T = 0JST c ). By definition, it is 
not symmetric under the reflection z <-> j8 — z. In the small z region, an effect from QQ + excited 
gluons can be seen, which is similar to the case of the standard hadronic correlation functions. In 
the intermediate z region, we find a single exponential behavior, which suggests the existence of 
a distinct peak of vanishing width in the spectral function. The position of the peak is nothing but 
the real part of the proper potential V(R,T). We make a single exponential fit of W\ for each R 
in the interval (15a T < z < 30a T ). The result is plotted as a function of R in the right panel of 
Fig.2 with filled squares. Error bars estimated by a % fit for several slightly shifted or contracted 
fitting regions reflect both the statistical and systematic uncertainties. The data can be fitted well 
by a Coulomb + linear form with an effective string tension y/o T=aiiTc = 367(5)MeV as indicated 
by the gray line. This is smaller than the known value at T = 0, ^jo T=a ~ 430MeV possibly due 
to the thermal fluctuations of the confining string. We emphasize here that (i) the thermal Wilson 
loop is directly linked to the Schrodinger equation as we have discussed and (ii) one can obtain not 
only the real part but also the imaginary part of the proper potential from the thermal Wilson loop 
by using e.g. MEM. Analysis of the present lattice data in the whole range of z by using MEM is 
currently under way. 

Shown in the left panel of Fig.3 is the thermal Wilson loop in the high temperature decon- 
finement phase (T = 2.33T C ). In this case, we find three characteristic features: For zja z <C N T , 
some evidence for excited states can be seen. For z/a t ~ N t /2, there appears an approximate single 
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Figure 3: (Left) The thermal Wilson loop W\(z : R) as a function of T at T = 2.33T C . The solid lines in the 
left panel are the results of the triple exponential fit in the intermediate region of T. (Right) Proper potential 
reconstructed from the triple exponential fit of the thermal Wilson loop. Color-singlet potential obtained 
from the Polyakov-line correlation in the color-singlet channel is also shown for comparison. 



exponential behavior which indicates a well-defined peak in the spectral function. For zja t ~ N t , 
the thermal Wilson loop increases as z increases, which indicates that the spectral function pi(cb,R) 
has some structure for negative ft) although the spectral strength is extremely small. This may be 
interpreted as the effect of thermal gluons with energy %; they can compensate a deficit of the 
external energy ft) to match the energy conservation, 2niQ = ft) + £ t h or equivalently ft) = — e±, so 
that a QQ can appear in the intermediate state. A triple exponential fit is deployed, designed to give 
stable results in the presence of non single-exponential behavior coming from the end regions of 
the z interval. The fit results of W\ for each R in the interval (10a T < z < 20a T ) leads to the plot 
in the right panel of Fig.3. It shows a Coulombic behavior at short distances, while the potential 
seems to be Debye screened at long distances. For comparison, we measure the color-singlet free 
energy on the same lattice from the Polyakov line correlations in the color-singlet channel under 
Coulomb gauge fixing: F^ l \R,T) = — rin(Trfl(x)n t (y)). The results are shown by the open 
circles. Although Re V(R,T) and F ( - 1 \R,T) have no direct theoretical connection, they coincide 
within the error margins devoid of any adjustments. An analysis of the present lattice data over the 
whole range of z by using MEM is currently under way to extract the complete spectral structure of 
Pi. This will enable us to extract lmV(R,T) and also to judge the validity of the potential picture 
at high T. 

4. Summary and concluding remarks 

We proposed a non-perturbative and gauge invariant approach to connect the Schodinger equa- 
tion description of a heavy QQ pair in terms of a static potential with the spectral function of the 
thermal Wilson loop obtained from Lattice QCD. It was shown that if the spectral structure is 
well defined for a set of temperatures (T) and separation distances (/?), the peak position and 
width correspond to the real and imaginary part of the proper potential V(R,T) respectively. A 
first determination of the real part of the proper potential from quenched lattice QCD simulations 
(Ng xN r = 20 2 x (96 - 32)) was presented for 0.787; and 2.33r c . It showed that although no 
apparent connection between the real part of the proper potential and the color-singlet free-energy 
potential exist, their values above T c coincide within the error bars. We are currently analyzing the 
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thermal Wilson loop data for 0.78 < T /T c < 2.33, using MEM to determine the full spectral struc- 
ture in an attempt to reconstruct not only the real but also the imaginary part, especially above the 
phase transition. In the case of infinitely heavy quarks, the only contribution to such an imaginary 
part at high T comes from the scattering of light medium particles, i.e. Landau damping [^|, |To| ]. 
The non-perturbative determination of Re V(R, T) and Im V(R, T) for 1 < T /T c < 2 is particularly 
important in relation to the fate of charmonium above T c . 

Two directions for future research are thus in order: (i) we need to include the effects of light 
fermions in the medium by utilizing full QCD Nj ■ = 2 + 1 configurations, and (ii) the approach 
has to be extended to finite QQ masses. As for (ii), a similar strategy as for T = 0, where the 
spatial fluctuations of the heavy quarks are incorporated by appropriate insertion of handles in the 
temporal Wilson lines [|12|] might be a first starting point. 

We thank Yuu Maezawa and the members of the WHOT-QCD Collaboration for useful dis- 
cussions. This research was supported in part by the Grant-in- Aid of MEXT (Nos. 18540253) and 
by Grant-in- Aid for Scientific Research on Innovative Areas (No. 2004: 20105003). 
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